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We consider in this paper a simple oscillating Quintom model of dark energy which has two free 
parameters and an equation of state oscillating and crossing -1. For low redshifts the equation of 
state of this model resembles itself similar to the linearly parameterized dark energy, however differ 
substantially at large redshifts. We fit our model to the observational data separately from the 
new high redshift supernova observations from the HST/GOODS program and previous supernova, 
CMB and galaxy clustering. Our results show that because of the oscillating feature of our model 
the constraints from observations at large redshifts such as CMB become less stringent. 



There are strong evidences that the Universe is spa- 
tially flat and accelerating at the present time 
In the context of Friedmann-Robertson- Walker cosmol- 
ogy, this acceleration is attributed to the domination of 
a component, dubbed dark energy. The simplest candi- 
date for the dark energy seems to be a remnant small 
cosmological constant or (true or false) vacuum energy, 
but it suffers from the difhculties associated with the fine 
tuning and coincidence problem, so many physicists are 
attracted with the idea that the acceleration is driven by 
dynamical scalar fields, such as Quintessence^!^. 

Despite the current theoretical ambiguity for the na- 
ture of dark energy, the prosperous observational data 
{e.g. supernova, CMB and large scale structure data and 
so on ) have opened a robust window for testing the re- 
cent and even early behavior of dark energy using some 
parameterizations for its equation of state. The recent 
fits to current supcrnova(SN) la data in the literature 
find that even though the behavior of dark energy is to 
great extent in consistency with a cosmological constant, 
an evolving dark energy with the equation of state w 
larger than -1 in the recent past but less than -1 today 
is mildly favored]^ Q- If such a result holds on with the 
accumulation of observational data, this would be a great 
challenge to current cosmology. Firstly, the cosmological 
constant as a candidate for dark energy will be excluded 
and dark energy must be dynamical. Secondly, the simple 
dynamical dark energy models considered vastly in the 
literature like the quintessence 0, 13 or the phantom*^ 
can not be satisfied either. 

In the quintessence model, the energy density and the 
pressure for the quintessence field are 

p^^Q^ + V{Q) , P=^Q'-V{Q) . (1) 

So, its equation of state w = p/p is in the range —1 < 
w < 1 for V{Q) > 0. However, for the phantom which 
has the opposite sign of the kinetic term compared with 
the quintessence in the Lagrangian, 

C^-^d^Qd'-Q-ViQ) , (2) 

the equation of state w = {—■^Q'^ — V)/{—^Q'^ + V) is 
located in the range of uj < —1. Neither the quintessence 



nor the phantom alone can fulfill the transition from w > 
— 1 to w < — 1 and vice versa. 

In Ref. we have proposed a new scenario of the dark 
energy, dubbed Quintom. A simple Quintom model con- 
sists of two scalar fields, one being the quintessence with 
the other being the phantom field. This type of Quin- 
tom model will provide a scenario where at early time 
the quintessence dominates with w > —1 and lately the 
phantom dominates with w less than —1, satisfying cur- 
rent observations A detailed study on the cosmologi- 
cal evolution of this class of Quintom model is performed 
in Ref. [Toj . The Quintom models are different from the 
quintessence or phantom in the determination of the evo- 
lution and fate of the universe. Generically speaking, the 
phantom model has to be more fine tuned in the early 
epochs to serve as dark energy today, since its energy den- 
sity increases with expansion of the universe. Meanwhile 
the Quintom model can also preserve the tracking behav- 
ior of q uintessence, where less fine tuning is needed. In 
Ref. nil we have studied a class of Quintom models with 
an oscillating equation of state and found that oscillat- 
ing Quintom can unify the early inflation and current 
acceleration of the universe, leading to oscillations of the 
Hubble constant and a recurring universe. Our oscil- 
lating Quintom would not lead to a big crunch nor big 
rip. The scale factor keeps increasing from one period 
to another and leads naturally to a highly flat universe. 
The universe in this model recurs itself and we are only 
staying among one of the epochs, in which sense the co- 
incidence problem is reconciled. In this paper we study 
the current observational constraints on the model of os- 
cillating Quintom. Speciflcally we consider a variation of 
the model in Ref. 1 1] with an equation of state given by 

w{z) — wq + wi sinz . (3) 

The model above has two parameters wq and wi and for 
low redshifts it reduces to 

w{z) = Wo + wiz , (4) 

which is commonly used in the literature for the linear 
parameterization of the dark energy. At large redshifts 
the w of our model oscillates and will differ substantially 
from the model with linear parameterization, which is 
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the key as we will show below to evade the constraints 
from the high redshifts observational data. 

We use the "ffold" set of 157 SN la data published by 
Riess et al. in "l] in fitting to supernova. We assume a 
flat space and use a uniform prior on the matter density 
fraction of the universe: 0.2 < ri,„ < 0.4. We constrain 
the Hubble parameter to be uniformly in 3a HST region: 
0.51 <h< 0.93. 

In a flat Universe, the luminosity distance can be ex- 
pressed as dL{z) = (1 + z)T{z)c/Ho, where r(z) = 
/o dz'/E{z') and 



E{z) = [an(l + zf + (1 - 



1/2 



Hiz)/Ho 



(5) 

where X{z) = px{z)/ px{0) and the label of "X" stands 
for dark energy. As an illustrative effect we use sim- 
ilar fitting methods by Refs.J^ For the Cosmic 
Microwave Background(CMB) fitting we use here only 

the CMB shift parameter i? = nll^T{zcMB) = 

1.716 ± 0.062 as in Ref.[l3 where ^cmb = 1089. The 
large scale structure(LSS) information we use is the lin- 
ear growth rate f{z2df ) = 0.51 ± 0.11 measured by the 
2dF survey (2dFGRS) [Hill, where Z2df = 0.15 and 
/ = {dlnD/d\na), D = 6k{z) / Sk{0) is the linear growth 
factor and the density perturbation equation of Sk evo- 
lutes as 6k + 2{a/a)5k — 4:TrGpMSk — 0. 




FIG. 1: SN la, CMB and LSS constraints on the oscillating 
dark energy model. The best fit values are shown in the cen- 
ters of each panel. The grey and light grey areas show the la 
and 2a confidence regions respectively. 
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FIG. 2; 1-dimensional likelihood of the variation amplitude 
of wi using SN la, CMB and LSS, shown with the la and 2a 
limits. 




FIG. 3: An illustration of the equation of state with redshift 
for best fit value of Figure □ using SN la + CMB + LSS, 
where {wo,wi) = (—1.33,1.47). 



In Figure U] we delineate SN la, CMB and LSS con- 
straints on the oscillating dark energy model. We can find 
a large parameter space is allowed by current supernova 
data while additional constraints from CMB R parameter 
and linear growth rate of 2df have restricted considerably 
the parameters within a lower region. We plot the cor- 
responding one dimensional likelihood of wi in Figure |2 
An oscillating Quintom model seems to be mildly favored 
at la. The best fit value for the SN la + CMB + LSS 
of our oscillating model is (wo,wi) = (—1.33,1.47). In 
Figure|31we show the corresponding w{z). As R contains 
the contribution of px on all relevant scales as an integra- 
tion, wi is somewhat constrained, this has also reflected 
the weak correlation of our parameterization. When LSS 
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FIG. 4: A comparison of linear growth factor and density 
between oscillating dark energy w — —1.33 + 1.47sinz and 

that with a constant equation of state w = -1.33. piQ q. a comparison of the maximum luminosity distance 

difference up to z = 1. for oscillating dark energy and the 
linear parameterized dark energy model. The contour stands 
for \AdL/dL\- 
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FIG. 5: SN la and CMB constraints on the linear parameter- 
ized dark energy model. Left : SN la only; Right: SN la + 
CMB constraints. 



constraint is also taken into account we find although the 
shape of the parameter space is not changed, the best fit 
parameter has a lower wi compared with SN la and CMB 
constraints. In the upper panel of Figure 0] we find that 
the linear growth factor deviates reasonably from a con- 
stant Wo case in low redshift. We should point out in 
our simple two parameter expansion the period of oscil- 
lation is fixed, which is due to the simplification of our 



parameterization. 

Before concluding we compare our model in (3) with 
the linearly parameterized dark energy model in (4). In 
Figure IHl we show the SN la and CMB constraints on lin- 
ear parameterized dark energy model in (4). Comparing 
Figure inland Figure ^ one can see the similarity for the 
linear and oscillating models when considering only the 
supernova data. This is understandable since our model 
in (3) for low redshifts reduces to the one in (4), how- 
ever when including the high redshift observational data 
such as CMB the difference is significant. To understand 
these firstly we show in Figure the maximum difference 
of Adi/di of the two dark energy models up to z < 1., 
where only 9 SN la of the gold sample lie outside the 
range, the maximum difference is no more than 2%. In 
the calculations we have set fim = 0.26 and h — 0.735 
as an example. The maximum difference is considerably 
small and the parameter spaces are expected to be similar 
for the two models under supernova constraints. Mean- 
while as from Figure[5lwe also expect wq ^ —0.5 for the 
oscillating case, in the early epoch as 

the contribution of wi term to px would be cancelled 
due to averaging over the oscillating periods, as can be 
seen from Figure 01 In Figure^ we have shown the ra- 
tio of density with w{z) — wq + wi sinz and w{z) = wq, 
assuming wq — —1.33 and Wi = 1.47 and the same back- 
ground parameters as Figure O We can see the ratio is a 
constant in the early epochs. Dark energy would hence 
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be negligible on large redshifts compared with matter en- 
ergy density and the R parameter of CMB would not give 
very strong constraints on the oscillating dark energy. 
When CMB is considered in the linear parameterization 
of the equation of state the upper limit shrinks promptly 
to near zero. This is because that px{z)/ px{0) would 
blow up for a sizable positive wi on large redshifts and 
hence CMB has restricted the model parameter wi much 
stronger than supernova. 

In summary we in this paper have considered a sim- 
ple two-parameter model of oscillating Quintom. The 
oscillating equation of state resemble itself similar to the 
linear parameterized dark energy in low redshift, mean- 
while the oscillating feature of the model make it pos- 
sible to evade the constraints from CMB observations. 
Our results indicate that a dynamical model of dark en- 



ergy such as the oscillating Quintom proposed in this 
paper is mildly favored under current observations. Our 
analysis here can be generalized into different Quintom 
model with an oscillating equation of state as a function 
of scale factor a for example ^/{a) — wq + w i sin (l — a) 
which becomes the one proposed by Refs.j^, [la for 
small (1 — a). 

Acknowledgements: We thank Alexei Starobinsky, 
Constantinos Skordis, Yun Wang, Jun'ichi Yokoyama 
and Zhong-Hong Zhu for helpful discussions. We thank 
Xiulian Wang, Xiaojun Bi, Peihong Gu, Hong Li and 
Yunsong Piao for useful conversations. This work is sup- 
ported in part by National Natural Science Foundation 
of China under Grant Nos. 90303004 and 19925523 and 
by Ministry of Science and Technology of China under 
Grant No. NKBRSF G19990754. 



[1] A. G. Riess et al, Astrophys. J. 607, 665 (2004). 

[2] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 

[3] U. Seljak et al., astro-ph/0407372 

[4] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 
(1988); P. J. E. Peebles and B. Ratra, Astrophys. ,]. 325, 
L17 (1988). 

[5] C. Wetterich, Nucl. Phys. B 302, 668 (1988); C. Wet- 
terich, Astron. Astrophys. 301, 321 (1995). 

[6] e. g. S. Nesseris and L. Perivolaropoulos, Phys. Rev. 
D 70, 043531 (2004); R. A. Daly and S. G. Djor- 
govski, astro-ph/0403664 U. Alam, V. Sahni, and 
A. A. Starobinsky, astro-ph/0403687 D. Huterer and 
A. Cooray, astro-ph/0404062 S. W. AUen, R. W. 
Schmidt, H. Ebeling, A. C. Fabian, and L. van Spey- 
broeck, astro-ph/0405340 T. R. Choudhury and T. 
Padmanabhan, astro-ph/0311622 S. Hannestad and 
E. MortseU, astro-ph/0407259 P. S. Corasaniti, M. 
Kunz, D. Parkinson, E. J. Copeland and B. A. Bas- 
sett, astro-ph/0406608 D. A. Dicus and W. W. Repko, 
astro-ph/0407094 B. A. Bassett, P. S. Corasaniti, and 
M. Kunz, astro-ph/0407364 _ 

[7] B. Feng, X. Wang and X. Zhang, astro-ph/0404224j 

[8] R. R. Caldwell, Phys. Lett. B 545, 23 (2002). 

[9] For relevant studies see e. g. S. Nojiri and S. D. 
Odintsov, Ph ys. Lett. B 59 9, 137 (2004); Y. H. Wei 
and Y. Tian, |gr-qc/0405038l C. Csaki and N. Kaloper, 



astro-p h/04095961 Y. Wang, J. M. Kratochvil, A. 
Linda, M. Shmakova, astro-ph/0409264 T. Brunier, 
V. K. Onemh, R. P. Woodard, gr-qc/0408080 B. 
Mclnnes, JHEP 0410, 018 (2004); G. Allemandi, A. 
Borowiec and M. Francaviglia, Phys. Rev. D 70, 103503 
(2004); E. Babichev, V. Dokuchaev and Y Eroshenko, 
lastro^ h/0407190 A. Vikman, astro-ph/0407107 Z. H. 
Zhu,"astro-ph/0404201 V. K. Onemh and R. P. Woodard 
gr-qc/0406098 S. D. H. Hsu, A. Jenkins and M. B. Wise, 
Phys. Lett. B 597, 270 (2004); Gang Chen and Bharat 
Ratra, Astrophys. J. 612, LI (2004); T. Padmanabhan, 
|astro-ph/0411044| 

[10] Z. Guo , T. Piao, X. Zhang and Y. Zhang, 
|astro-ph /0410654 

[11] B. Feng, M. Li, Y.Piao and X. Zhang, astro-ph/0407432j 

[12] R. A. Knop et al, Astrophys. J. 598, 102 (2003). 

[13] Y. Wang and M. Tegmark, Phys. Rev. Lett. 92, 241302 
(2004). 

[14] J. R. Bond, G. Efstathiou and M. Tegmark, MNRAS 

291, L33 (1997). 
[15] E. Hawkins et al, MNRAS 346, 78 (2003); L. Verde et 

al, MNRAS 335, 432 (2002). 
[16] M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 

213 (2001). 

[17] E. Linder, Phys. Rev. Lett. 90, 091301 (2003). 



